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We address the reconstruction of the full photon distribution of multimode fields generated 
by seeded parametric down-conversion (PDC). Our scheme is based on on/off avalanche pho- 
todetection assisted by maximum-likelihood (MaxLik) estimation and does not involve photon 
counting. We present a novel constrained MaxLik method that incorporates the request of 
finite energy to improve the rate of convergence and, in turn, the overall accuracy of the 
reconstruction. 



Introduction 



The reconstruction of photon statistics of quantum optical fields is of the utmost 
relevance for several applications, ranging from quantum information (y|) to the 
foundations of quantum mechanics H) and quantum optics (0, 0, [3, [3, 0)- Despite 
this fact, the realization of photodetectors well suited for this purpose still repre- 
sents an experimental challenge. The few existing examples {3, [3. [73. 11 , li. 



show severe limitations. On the other hand, reconstruction schemes based on quan- 
tum tomography ( 1^, 73, 73) require phase-matching with a suitable local oscillator 



and do not represent a technique suited for a diffuse use. This situation prompted 
various theoretical studies (fTj . [73, [73, (13) addressed to achieve the reconstruc- 
tion of the (diagonal) elements of the density matrix exploiting the information 
achievable with realistic detectors. 
In a recent series of papers {21, H, 13, 3, H)> we have demonstrated how a very 



satisfactory reconstruction of the statistics of mono-partite and bi-partite quantum 
optical states may be obtained using the simplest kind of detectors, namely on/off 
detectors C3, \2(h operating in the Geiger mode, whose outcomes are either "off" 
(no photons detected) or "on" , i.e., a "click" , indicating the detection of one or more 
photons. Our method recovers the full photon statistics using maximum likelihood 
(MaxLik) reconstruction on on/off data obtained using variable detection efficiency 
(by inserting calibrated neutral filters). 

In this paper we present a modified version of our method, which improves both 
convergence and accuracy upon incorporating the obvious a priori constraint of 
finite signal energy. In particular, we address the reconstruction of the full photon 
distribution of multimode fields generated by seeded parametric down-conversion. 
The novel method allows to overcome the increased complexity of the reconstruc- 
tion problem and represents an important step in view of a widespread application 
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of MaxLik reconstruction. 

The paper is structured as follows: in the next Section we describe the constrained 
MaxLik algorithm in some details, whereas in Section[3]we describe the experimen- 
tal apparatus and illustrate the application of the method to the reconstruction of 
the photon distribution of multimode fields generated by seeded parametric down- 
conversion. Finally, Section [4] closes the paper with some concluding remarks. 



2. Constrained MaxLik algorithm 

The probability po(r]) that a photodetector with quantum efficiency r] does not 
click, when an input quantum state g = Y^ n m Qnm\n)(m\ impinges on it, reads as 
follows: 

p (r,)=J2^-v) n Qn, (1) 

n 

where g n = g nn is the n-th entry of the photon distribution of the input state. 
Now, if we consider a set of N detectors with different quantum efficiencies r) v , 
v = 1, . . . , N, then we can write the "off" probabilities as: 

P„ = p {r) u ) = y y A vn g n , (2) 

n 

where A vn = (1 — r/ u ) n . Looking at Eq. ([2]) as a statistical model for the parameters 
g n , we can solve it by the MaxLik estimation. We proceed as follows: first of all, 
we assume there exists a value h such that g n is negligible for n > h; we assign the 
loglikelihood function (with normalized P v ), that is the global probability of the 
sample: 

where f v = M u /N x is the experimental frequency of "off" events, M u being the 
number of "off" events for a fixed quantum efficiency rj u and N x the total events 
amount. The MaxLik estimated g n values are the ones maximizing L. Since the 
model is linear and the unknowns g n are positive, we can find the expectation- 
maximization solution of the MaxLik problem by means of an iterative procedure 



as described in \12, \M, \M, \2l, \M), i.e.: 



n W A f 

where gn is the value of g n evaluated at the h-th iteration, and = ^ n A un g^"\ 
The algorithm (JU), known to converge unbiasedly to the MaxLik solution, provides 
a solution once the initial distribution g^ is chosen. On the other hand, the initial 
distribution slightly affects only the convergence rate and not the precision at 
convergence (jid ). 

As a matter of fact, the solution obtained above corresponds to the best photon 
distribution fitting the experimental data, i.e., the measured "off" probabilities. 
However, it is possible that different photon distributions fit the same experimental 
data, giving rise to a family of suitable distributions. In these cases it could happen 
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that the MaxLik solution, even if in good agreement with the available experimental 
knowledge, may be different from the actual (unknown) one. Indeed, this is the 
case of multimode fields when the number of modes grows. In order to overcome 
this limitation, we developed a modified version of the MaxLik algorithm, which 
incorporated the constraint of finite energy for the incoming signal, i.e., the quantity 
^2n n 8n- In practice, we maximize the loglikelihood ([3j) with a constraint on the 
energy. Now the function to be maximized with respect to g n is: 



nQn, (5) 



L being given in ([3]) and (3 being a Lagrange multiplier. The equations ^ 
lead to: 

£ : — /e,pvm p- = 1 ' ^ 



and, then, by multiplying both the sides of Eq. © by g n , we get a map Tg n = g 
whose fixed point can be obtained by the following iterative solution: 



(fc+l) = Un 

2_jm fm v 



o {h) ^ A f 



(7) 



The parameter (3 can be tuned in order to control the energy of the reconstructed 
state and improve both the convergence rate and the overall accuracy. Of course, 
if we take (3 = 0, then Eq.s ([7]) and (JH) become the same. 

Indeed, in order to use Eq. ([7]) we need to know the input state energy which, 
in general, cannot be directly accessible from experimental data. However, in cases 
when a model of the photon distribution of the input state is available, we can 
estimate indirectly this energy by a simple fit of the "off" probabilities. 

In the following, we consider the multimode field obtained by seeded Parametric 
Down Conversion (PDC). In this case, the photon distribution is expected to have 
the form: 



™ (1 + AW" +M V 1 + AW " V Wthtl+JVth 

where g n = g n (N t h,a, M) and L®(z) are the generalized Laguerre polynomials. 
Eq. (J8j) represents the convolution of M thermal states (M spatial modes), all with 
the same average number of thermal photons -/V tn > except for one, displaced 
by an amount a. This model will be justified by the experimental setup described 
in Section [3l From Eq.s ([1]) and (JSj) we can calculate the "off" probability po = 
p (N th ,a,M,rj), that is: 

00 1 / I p \ 

Po = £(1 - I)" 9n(N th , a, M) = exp (-jf^) , (9) 



r] being the quantum efficiency of the on/off photodetector. Notice that, from 
Eq. ([9]) , one can obtain the relevant cases for a Poissonian input photon distribution 
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(-^th —* 0) and for a multithermal one (a — * 0). Thanks to Eq. ([9]), one can evaluate 
the input state energy iV t h + |a| 2 , and choose a suitable value for j3. 

Before the end of this section, it is worth pointing out that if each spatial mode 
consists of M' temporal modes, then the input photon distribution and the "off" 
probability are still given by Eq.s (jHj) and ([9]), respectively, but with M->Mx M' 
in place of M. 



3. Experimental test 

In order to test the reliability of the algorithm reported in Eq. ([7]), we applied it 
to the reconstruction of a stimulated type-I PDC branch, at different stimulation 
regimes. In our experiment, whose setup is shown in Fig. [H a CW Argon laser 
(Apump = 351.1 nm) pumps a 5 x 5 x 5 mm type-I BBO crystal, generating PDC. 
Together with the pump beam, a CW Nd:Yag laser (A see d = 1064 nm) is injected 
into the crystal in the proper way to generate stimulated PDC, and we look at the 
emission in the k s ti mu i direction (A s timul = 524 nm). Different values of quantum 
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Figure 1. Experimental setup: the stimulated emission at A st ; mu i = 524 nm is addressed to the NF and 
then collected by the APD. The number and temporal width of the acquisition windows is set by the pulse 
generator used for the detector's gating. 

efficiency have been obtained by inserting properly calibrated Schott neutral filters 
(NF), starting from 7y max = 28.4%; after them, on the optical path of the stim- 
ulated branch we put an anti-infrared (IR) filter (to cut off the noise due to the 
Nd:Yag laser dispersion), a variable pinhole (to control the number M of spatial 
propagation modes collected) and a fiber coupler connected by a multimode fiber 
with the detector (avalanche photodiode, Perkin Elmer SPCM-AQR-15). We set 
the pulse generator in order to open in the APD 2 10 5 detection windows per sec- 
ond, each one of 20 ns; the pinhole diameter is regulated in order to collect only few 
spatial modes (more precisely M = 7), of whom only one stimulated. Moreover, 
each spatial mode consists of many temporal modes: the total number of modes 
can then be estimated M = 710 5 . It is worth to mention that, when the number 
of modes exceeds few tens, the dependence on this parameter is rather small and 
a rough estimate of the order of magnitude suffices. 

We have performed three separate data collections, each one corresponding to 
a different stimulation regime: by indicating with x the percentage of stimulated 
emission on the whole PDC amount collected, our acquisitions were respectively 
characterized by x = 51.4%, x = 78.1% and x = 90.7%. The evaluation of the 
background photons have been performed through an acquisition step without PDC 




Figure 2. On the left: fo non-click frequencies (gray disks) given by the stimulated PDC with different 
stimulation regimes as functions of the quantum efficiency r\. The black disks are the "off" probabilities 
obtained by means of the MaxLik reconstructed photon distribution; the solid line corresponds to Eq. J9} 
with jV t h = (1 — x)N avc , ct | 2 = xiV a ve and M = 7 10 5 . In each plot we report also the average number of 
photons (./V ave ) obtained by the fit of the experimental fo, the percentage of stimulated emission x and the 
X 2 of the the MaxLik fit. On the right: MaxLik reconstructed photon distribution (gray bars) and photon 
distribution given by Eq. JS} with the same values of the parameters given in the respective left plots. In 
each plot we report also the fidelity F between the two photon distributions. Note the different ranges of 
n. 



The obtained results are shown in Figure EJ for the reconstruction we used the 
MaxLik estimation with constraint on the energy, as described in the previous 
section. The plots on the left show the /o non-click frequencies given by the stim- 
ulated PDC with different stimulation regimes vs. the quantum efficiency 77, and 
the fit obtained by means of the MaxLik estimation and Eq. ([9]). The x 2 quantity 
reported has been defined as the sum of the square differences between the "off" 
probabilities given by the reconstructed photon statistics and the measured fo. To 
quantify the similarity between the two photon distributions appearing in the plots 
on the right, instead, we used the fidelity formula: 

f = E^ ml) ^ dmt) , do) 

n 

where g^ L ^ is the MaxLik reconstructed photon distribution and g ( h >MT ^ is the 
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one obtained from Eq. ([8]) . In Figure [3] we consider the same scenario giving the 
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Figure 3. Same plots as the bottom ones in Figure [2] obtained without the energy constraint on the 
algorithm (notice that here AT ave = 17.55 whereas previously it was 7V ave = 18.25): even if the MaxLik fit 
of the experimental /o set is good, the fidelity between the MaxLik reconstructed photon distribution and 
the one given by Eq. |J5} is quite low. 

bottom plots of Figure [2j but now we try to perform the reconstruction without 
any constraint on the energy (/? = 0). We can see that, even if the MaxLik fit of 
the fo frequencies (left side) is good, the fidelity between the MaxLik reconstructed 
photon distribution and the one given by Eq. ([8]) is quite low (right side): a result 
that confirms the advantage of using the constrained MaxLik method. 



4. Concluding remarks 

In this paper we have shown how an important improvement on the convergence 
of the photon statistics reconstruction code, based on MaxLik estimation applied 
to on/off detection data, can be achieved by increasing the number of Lagrange 
multipliers when some "a priori" knowledge of the state is available. In particular we 
have addressed the reconstruction of the full photon distribution of multimode fields 
generated by seeded parametric down-conversion, demonstrating the advantages of 
the constrained MaxLik method. This achievement represents an important step 
in view of widespread application of this scheme. 



Acknowledgements 

This work has been partially supported by the CNR-INFM convention, by Regione 
Piemonte E14 contract and by 07-02-91581-ASP. 



References 



(1) Bouwmeester, D.; Ekert, A.K.; Zeilinger, A. The Physics of Quantum Information: Quantum Cryp- 
tography, Quantum Teleportation, Quantum Computation, Springer, New York, 2000. 

(2) Genovese, M. Physics Reports, 2005, 413/6. 

(3) Vogel, W.; Welsch, D.G. Quantum Optics - An Introduction, Wiley-vch, Berlin, 2001. 

(4) Schleich, W.P. Quantum Optics in Phase Space, Wiley- Vch, Berlin, 2001. 

(5) Perina, J.; Hradil, Z.; Jurco, B. Quantum Optics and Fundamental Physics, Kluwer, Dordrecht, 1994. 

(6) Leonhardt, U. Measuring the Quantum State of Light , Cambridge Univ. Press, Cambridge, 1997. 

(7) Mandel, L.; Wolf, E. Optical Coherence and Quantum Optics, Cambridge Univ. Press, Cambridge, 
1995. 

(8) Zambra, C; Bondani, M., Rev. Sci. Instrum. 2004, 75, 2762-2765. 

(9) Kim, J.; Takeuchi, S.; Yamamoto, Y. Appl. Phys. Lett. 1999, 74, 902-904. 

(10) Peacock, A.; Verhoeve, P.; Rando, N.; van Dordrecht, A.; Taylor, B.C.; Erd, C; Perryman, M.A.C.; 
Venn, R.; Howlett, J.; D. J. Goldie, D.J.; Lumley, J.; Wallis, M. Nature 1996, 381, 135-137. 

(11) Zappa, F.; Lacaita, A.L.; Cova, S.D.; Lovati, P. Opt. Eng., 1996, 35, 938-945. 

(12) Achilles, D.; Silberhorn, C; Liwa, C; Banaszek, K.; Walmsley, LA. Opt. Lett. 2003, 28, 2387-2389. 



February 15, 2008 



9:57 



Journal of Modern Optics 



stimul'papcrl 



REFERENCES 



7 



(13) Di Giuseppe, G.; Sergienko, A.V.; Saleh, B.E.A.; Teich, M.C. Quantum Information and Computa- 
tion, Proc. SPIE, 2003 5105, 39-50. 

(14) Munroe, M.; Boggavarapu, D.; Anderson, M.E.; Raymer, M.G. Phys. Rev. A , 1995, 52, 924-927. 

(15) Zhang, Y.;, Kasai, K.; Watanabe, M. Opt. Lett, 2002, 27, 1244-1246. 

(16) Raymer, M.; Beck, M. Quantum States Estimation, Lect. Not. Phys., 649, 235-295, Springer, Berlin- 
Heidelberg, 2004. 

(17) Mogilevtsev, D. Opt. Coram, 1998, 156, 307-310. 

(18) Mogilevtsev, D. Acta Phys.Slov., 1999, 49, 743-748. 

(19) Rossi, A.R.; Olivares, S.; Paris, M.G. A. Phys. Rev. A, 2004, 70, 055801. 

(20) Rossi, A.R.; Paris, M.G. A. E. Phys. Jour. D, 2005, 32, 223-226. 

(21) Zambra, G.; Andrconi, A.; Bondani, M.; Gramcgna, M.; Genovese, M.; Brida, G.; Rossi, A.; Paris, 
M.G. A. Phys. Rev. Lett., 2005, 95, 063602/1-4. 

(22) Gramegna, M.; Genovese, M.; Brida, G.; Bondani, M.; Zambra, G.; Andreoni, A.; Rossi, A.R.; Paris, 
M.G. A. Laser Physics, 2006, 16, 385-392; 

(23) Brida, G.; Genovese, M.; Gramegna, M.; Paris, M.G. A.; Predazzi, E.; Cagliero, E. Open Systems & 
Information Dynamics, 2006, 13, 333-341. 

(24) Brida, G.; Genovese, M.; Piacentini, F.; Paris, M.G. A. Optics Letters, 2006, 31, 3508. 

(25) Brida, G.; Genovese, M.; Paris, M.G. A.; Piacentini, F.; Predazzi, E.; Vallauri, E. Optics and Spec- 
troscopy, 2007, 103, 95. 

(26) Banaszek, K. Phys. Rev. A, 1998 57, 5013. 

(27) Hradil, Z.; et al., Phys. Rev. Lett., 2006, 96, 230401. 

(28) Olivares, S.; et al., EPL, 2007, 80, 64002. 



